Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease.

The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b, the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.

Molecular Docking and Virtual Screening Based Prediction of Drugs for COVID-19.

AIMS: To predict potential drugs for COVID-19 by using molecular docking for virtual screening of drugs approved for other clinical applications. BACKGROUND: SARS-CoV-2 is the betacoronavirus responsible for the COVID-19 pandemic. It was listed as a potential global health threat by the WHO due to high mortality, high basic reproduction number, and lack of clinically approved drugs and vaccines. The genome of the virus responsible for COVID-19 has been sequenced. In addition, the three-dimensional structure of the main protease has been determined experimentally. OBJECTIVE: To identify potential drugs that can be repurposed for treatment of COVID-19 by using molecular docking based virtual screening of all approved drugs. METHODS: A list of drugs approved for clinical use was obtained from the SuperDRUG2 database. The structure of the target in the apo form, as well as structures of several target-ligand complexes, were obtained from RCSB PDB. The structure of SARS-CoV-2 Mpro determined from X-ray diffraction data was used as the target. Data regarding drugs in clinical trials for COVID-19 was obtained from clinicaltrials.org. Input for molecular docking based virtual screening was prepared by using Obabel and customized python, bash, and awk scripts. Molecular docking calculations were carried out with Vina and SMINA, and the docked conformations were analyzed and visualized with PLIP, Pymol, and Rasmol. RESULTS: Among the drugs that are being tested in clinical trials for COVID-19, Danoprevir and Darunavir were predicted to have the highest binding affinity for the Main protease (Mpro) target of SARS-CoV-2. Saquinavir and Beclabuvir were identified as the best novel candidates for COVID-19 therapy by using Virtual Screening of drugs approved for other clinical indications. CONCLUSION: Protease inhibitors approved for treatment of other viral diseases have the potential to be repurposed for treatment of COVID-19.